Symmetry preserving lattice collapse in tetragonal SrFe2- :r R<u a ;As2 (x = 0, 0.2) — 
a combined experimental and theoretical study 
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In a joint experimental and theoretical study, we investigate the isostructural collapse from the 
ambient pressure tetragonal phase to a collapsed tetragonal phase for non-superconducting metallic 
SrFe2As2 and SrFe1.gRuo.2As2. The crystallographic details have been studied using X-ray powder 
diffraction up to 20 GPa pressure in a diamond anvil cell. The structural phase transition occurs 
at 10 GPa and 9 GPa for SrFe2As2 and SrFe1.8Ruo.2As2, respectively. The changes in the unit 
cell dimensions are highly anisotropic with a continuous decrease of the c lattice parameter with 
pressure, while the a-axis length increases until the transition to a collapsed tetragonal phase and 
then continues to decrease. Across the phase transition, we observe a volume reduction of 5% and 
4% for SrFe2As2 and SrFe1.sRuo.2As2, respectively. We are able to discern that Ru substitution on 
the Fe-site acts like 'chemical pressure' to the system. Density-functional theory-based calculations 
of the electronic structure and electron localizability indicator are consistent with the experimental 
observations. Detailed analysis of the electronic structure in fc-space and real space reveals As 4p z 
interlayer bond formation as the driving force of the c/a collapse with a change in the As- As bond 
length of about 0.35A. 
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I. INTRODUCTION 



The surprising discovery of superconductiv- 
ity in Fe-based pnictides and chalcogenides has 
steered and revived the interest in high-temperature 
superconductivity^ Most of the superconducting mem- 
bers discovered in the past 3 years can be grouped 
into five families named by their parent compounds: 
(i) 1111 family - i?OFeAs / AFeAsF (R = rare-earth, 
A=alkaline-earth metal); (ii) 122 family - AFc 2 As 2 (A 
= alkaline-earth or divalent rare-earth metal); (hi) 111 
family - ^4FeAs (A = alkali- metal) ; (iv) 11 family - 
FeSe(Te); (v) P22 family - PFc 2 As 2 (P = perovskite- 
oxide like Sr4Sc 2 0g). The basic common feature of 
these new parent compounds is the FeAs building block 
separated by spacer layers comprising of the above men- 
tioned alkali, alkaline-earth, rare-earth oxide/fluoride or 
a perovskite-oxideii~— Superconductivity is obtained by 
suitable doping of the parent compounds. For selected 
recent exhaustive overviews, we refer the readers to 
Refs.0 and0. The thickness of the spacer layers governs 
the extent of the quasi-two dimensional (2D) nature of 
the electronic structure. Using plasma frequencies as 
a tool to identify the 'effective dimensionality' among 
these systems, it has been shown that the 1111 systems 
are considerably more 2D than the 122 systems^ 
Interestingly, within these two families a correlation 
between the dimensionality and the superconducting 
transition temperature (T c ) can be construed, with the 
more 2D-like 1111 systems exhibiting a larger T c than 
the less anisotropic 122 systems^ Superconductivity in 
these systems upon the suppression of spin-density-wave 
(SDW) antifcrromagnctic order of Fc can be realized 
via hole-doping, electron-doping, isovalent substitution, 



or pressure £ut~— In contrast to doping, high pressure 
techniques provide a cleaner route to modify the elec- 
tronic structure without the added effects of chemical 
complexity^ Many of the high-pressure experiments 
were focused on investigating the suppression of the 
Fc-SDW order and the enhancement of T c . Recently, 
pressure-induced isostructural transitions from a param- 
agnetic tetragonal phase (T) to a collapsed tetragonal 
(cT) phase have been reported for several members of 
the 122 family: CaFe 2 As 2 , BaFe 2 As 2 , EuFe 2 As 2 , and 
EuCo 2 As 2 (Refs.Hl-ill)- The high electronic flexibility 
of compounds with the ThCr 2 Si 2 -type crystal structure 
can be related to the subtle interplay of covalent, ionic 
and metallic bonding contributions as has been outlined 
in several previous studies^r— Hoffmann and Zheng22 
also pointed out that in certain cases this might lead 
to structural phase transitions (see below). Indeed, 
these structural phase transitions of first and second 
order have been observed in many AT 2 P 2 compounds 30 
between two crystallographically isostructural modifica- 
tions with significantly different P-P distances along the 
tetragonal c axis. Common to them are the dramatic 
and highly anisotropic changes in unit cell dimensions 
(changes of the tetragonal c parameter by -10% and 
counteracting changes of a by +2%, leading to a collapse 
of the unit cell volume) depending on temperature^ 
pressure ; 30 ' 33 " — composition^ and chemical pressured 
Curiously, the order of the phase transition (first or 
second-order) was suggested to be dependent on the 
transition metal atom (LaFe 2 P 2 - first-order, LaCo 2 P 2 - 
second-order)^ - In contrast, recent experiments report 
the transition to be of second order for the arsenides 
EuX 2 As 2 (T = Fe, Co). More systematic studies are 
necessary to clarify the nature of these transitions, 
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as well as to understand the underlying physics that 
accompany the formation of the cT-phasc. 

In our joint theoretical and experimental study on 
SrFe2As2 and the isovalent Ru substituted sample 
SrFe1.8Ruo.2As2, we attempt to address a multitude of 
issues. First we analyze the possibility of a T — > cT 
phase transition in these systems as has been suggested 
for many other compounds with ThCr2Si2-type crys- 
tal structure. Previous theoretical studies&21 have pre- 
dicted a phase transition from the magnetically-ordered 
orthorhombic SDW phase to a tetragonal phase, but the- 
oretical studies on T — > cT phase are rather cursory. 
Secondly, we address the nature of the chemical bonding 
in the 122-arsenides and compare it to the well stud- 
ied 122-phosphides. Finally, we challenge the concept of 
chemical pressure by isovalent substitution of Fe by Ru 
and study the nature of the isostructural phase transi- 
tion and the bonding situation of the pnictide. For an 
isovalent substitution on the Fe site, the As layers re- 
main intact and unaffected from impurities arising from 
the substitution elements. 



II. METHODS 

Polycrystalline samples have been synthesized using 
solid-state reactions, similar to those described in Refs.@ 
and [l2l . Samples were obtained in the form of sintered 
pellets. X-ray diffraction measurements (XRD) were per- 
formed at the high-pressure beam-line ID09 of the ESRF 
up to 20 GPa at room temperature for SrFc2As2 and 
SrFe1.gRuo.2As2. For best possible hydrostatic condi- 
tions we used a membrane diamond anvil cell (DAC) with 
helium as pressure transmitting medium. The pressure 
was determined using the ruby fluorescence method. The 
measured powder rings were integrated using the pro- 
gram FIT2Di2£ After a background correction the lattice 
parameters were determined with the FullProi^ package. 

Density-functional theory (DFT) based band-structure 
calculations were performed using a full-potential all- 
electron local-orbital code FPL O 40 ' 41 within both local 
density approximation (LDA) and as well as generalized 
gradient approximation (GGA). Rclativistic effects were 
incorporated on a scalar-relativistic level. A well con- 
verged fc-mesh with 24 3 points in the full Brillouin zone 
was used. The crystal structures were optimized at dif- 
ferent levels to investigate or isolate effects that may de- 
pend sensitively to certain structural features. The full- 
relaxation of the unit cell involves optimizing the c/a 
ratios in addition to relaxing the z(As) coordinate. The 
electron localizability indicator /function (ELI/ELF) was 
evaluated according to Ref. [42] with an ELI/ELF mod- 
ule implemented within the FPLO program packaged 
The topology of ELI was analyzed using the program 
Basin— with consecutive integration of the electron den- 
sity in basins, which are bound by zero-flux surfaces in 
the ELI gradient field. This procedure, similar to the one 
proposed by Bader for the electron density— allows to 



assign an electron count for each basin. 

III. RESULTS AND DISCUSSION 
A. Symmetry preserving lattice collapse 

1. Experiment: SrFe2As2 

Collected in Fig. Q] are the measured tetragonal lat- 
tice parameters a and c, the axial ratio c/a and the unit 
cell volume as a function of pressure for the parent com- 
pound SrFe2As2 (the 10% Ru-doped SrFc1.gRuo.2As2 will 
be addressed later in this report). For SrFe2As2, the lat- 
tice parameters a and c show an initial decrease with 
pressure up to 5 GPa. Upon further increase in pres- 
sure, anomalous compression effects are observed with 
the lattice parameter a expanding rapidly up to around 
10 GPa, while the c lattice parameter continues to de- 
crease. At this juncture a strong decrease is witnessed 
in the c lattice parameter to smaller values (a transition 
from a T — >cT phase), after which a normal compres- 
sion behavior is observed up to 20 GPa. The c/a ratio as 
a function of pressure shows the onset of the structural 
phase transition around 10±1 GPa from a T-phase (with 
c/a « 2.85) to a cT-phase (with c/a w 2.65). Please re- 
fer to Table. 1 and 2 in the supplementary information 
for more details. 

Similarly, an abrupt change in the interlayer As-As 
distance was observed. However, this change of the 
interlayer distance can be acountcd for only partially 
by the change of the lattice parameter c (see Fig. [U 
lower panel). This distance decreases from the rather 
large value of 3.41 A (at ambient pressure) to a value of 
2.59 A (at the highest measured pressure of 19.6 GPa). 
Similar interatomic distances are found for example in 
skutterudite CoAs 3 (2.49/2.56 A)2£ or filled skutterudite 
LaFe 4 Asi2 (2.57/2.58 A) 41 These values are only slightly 
larger than the single-bond distance of 2.52 A occurring 
in a-As4£ The collapse of the lattice is in turn, also 
accompanied by a simultaneous pronounced decrease of 
the Sr-As distances ds r As and Fe-As distances dFoAs by 
about 12% and 4%, respectively, as well as the angles in 
the FeAs4 tetrahedron (see supplementary information). 
Electronic structure calculations have shown that even 
minute shifts in z (As) and thus the resulting Fe-As dis- 
tances have strong impact on the occupation of the Fe 
3d x 2_ y 2 orbitals, and therefore, the magnetic behavior 
in the ^4Fc2As2 phases . 21 ' 49 These effects should become 
even more magnified during the observed phase transi- 
tion. In contrast, under applied external pressures up to 
6 GPa, a pronounced robustness of the Fe-As bonds is 
recognized for BaFe2As2^ which crystallizes with a sig- 
nificantly larger c/a ratio of 3.35. However, the Ba-As 
distances significantly contract (about 3%). Neverthe- 
less, no indications of an onset of a lattice collapse were 
detected for this compound. 

The measured pressure-volume data were fit using a 
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FIG. 1. (Color Online) The measured tetragonal lattice 
parameters a and c, the axial ratio c/a and the unit cell 
volume as a function of applied pressure for SrFe2As2 and 
SrFe1.8Ruo.2As2. The measurements were performed at room 
temperature. The inset shows the Birch-Murnaghan equation 
of state fitting for the tetragonal (T) and collapsed (cT) phase 
for SrFe2As2. The lowest panel displays As- As distance for 
SrFe2As2 obtained from a full refinement of the z coordinate 
of the As atoms. To separate the influence due to the decrease 
in c, only, the As- As distance using the z coordinate of As at 
ambient pressure is also shown. 



Birch-Murnaghan equation of state for the two separate 
phases (shown as an inset in Fig. [1]). The obtained equi- 
librium volume Vq, bulk modulus Bq and its pressure 
derivative B' for the tetragonal phase (0 - 8.6 GPa) are 
V = 190.7 A 3 , B = 63.63 GPa and B' = 2.51, respec- 
tively, and for the collapsed tetragonal phase (11 - 19.5 
GPa) are V = 180.8 A 3 , B = 89.57 GPa and B' = 2.51, 




FIG. 2. (Color Online) Calculated a, c and c/a as a function 
of unit cell volume for SrFe2As2 within LDA. At each volume, 
only the c/a optimization has been performed. The c lattice 
parameter decreases throughout, while the a lattice parameter 
undergoes an anomalous expansion for some reduced volumes 
before it continues to decrease again. The vertical line denotes 
the experimental ambient condition volume. 

respectively. Our observation of the T — > cT-phase tran- 
sition in SrFe2As2 along with the anisotropic compress- 
ibility phenomena is akin to the reports on other 122 
arsenides AT 2 As 2 {A = Ca, Ba, Eu; T = Fe, Co)^Si 

2. Calculations: SrFc2As 2 

Previous work on the tetragonal collapse in other mem- 
bers of the AFe2As2 family (A = Ca, Ba, and Eu) have 
mostly been experimental studies*^— Analysis of the 
electronic structure of the Fe-As-based systems using 
band structure calculations has been quite successful in 
describing certain 'general trends' like the tetragonal to 
orthorhombic structural distortion, reduction in the mag- 
netic moment as a function of doping, etc., though DFT 
is less accurate in reproducing certain details like the 
Fe-As bond length.— In 2009, Yildirim reported on re- 
sults from first-principles calculations for CaFe2As2 and 
noted that the formation of the cT-phasc is controlled 
by the Fe-spin stated He suggested that pressure re- 
duces the Fe-spin moment which in turn weakens Fe-As 
bonding and strengthens As- As interaction and therefore 
causes the collapse of the lattice parameters. Here, also 
using first-principles calculations, we de-couple Fe mag- 
netism from the chemical bonding scenario and elucidate 
the bond formation as the key feature. To realize this, 
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FIG. 3. (Color Online) Calculated a, c and c/a as a function 
of unit cell volume for SrRu2As2 within LDA. At each volume, 
only the c/a optimization has been performed. The vertical 
line denotes the experimental ambient condition volume. 



we calculated the change in the ground state energy as a 
function of volume employing both LDA and GGA and 
using the tetragonal symmetry without invoking spin de- 
grees of freedom explicitly. At each volume two kinds of 
optimization were performed, one in which only the c/a 
ratio was optimized and another in which a simultane- 
ous optimization of the z(As) position was additionally 
carried out. Collected in Fig. [5] are the a and c lattice 
parameter and the c/a ratio as a function of volume ob- 
tained using the LDA with only c/a optimization (refer 
to Fig. 1. in the supplementary information for details 
regarding the GGA calculations and the effects of z(As) 
optimization). The calculations reproduce the experi- 
mental results including the transition from the ambient 
pressure T-phase to the cT-phase for SrFe2As 2 in accord 
with a dramatic jump in the axial ratio c/a from 2.85 
—> 2.65. Our calculations also reproduce the observed 
anomalous expansion of the lattice parameter a under 
pressure. This is the first theoretical observation of the 
anomalous expansion phenomena for the FcAs family of 
materials. By fitting the energy-volume (E— V) curves to 
two separate Birch-Murnaghan equations of state (EOS) 
fits (not shown here), we obtain a transition pressure 
of 11±1 GPa. The pressure (P) is obtained from the 
volume derivative of the EOS, which is inverted to get 
V{P). Equating the enthalpies E[V(P)] + PV(P) of the 
two phases gives the transition pressure. The calculated 
value is in good agreement with the experimental value 
of 10±1 GPa. Calculations using GGA (refer to Fig. 1. 



in supplementary information) with only c/a optimiza- 
tion give essentially the same trend for a, c and c/a as 
described above for LDA. 

Relaxing additionally the z(As) position worsens the 
quantitative description for both LDA and GGA. This 
behavior is not surprising because, as mentioned previ- 
ously, the Fe-As bond length is quite sensitive and DFT 
fails in accurately reproducing the experimental Fe-As 
bond length even for ambient conditions ! 52 ' 53 Another 
inadequacy of DFT that is worth mentioning here is the 
equilibrium volume. It is well known that LDA and GGA 
usually underestimate and overestimate, respectively, the 
equilibrium volume by just a few percent as compared to 
experiments. For SrFe2As2 both LDA and GGA underes- 
timate the equilibrium volume by 17% and 10%, respec- 
tively. This is quite unusual and at present there exists 
no conclusive reasoning for such a behavior. One pos- 
sible explanation that is gaining more acceptance is the 
presence of a nematic order, which is unaccessible using 
the present-day DFT toolsi 5 ^ 



3. Experiment: SrFe1.gRuo.2As2 

Similar to hydrostatic external pressure, studies on 
the effects of substitution on the transition metal Fe- 
site also show anisotropic changes in the lattice param- 
eters, with a significant contraction of the c-axis length 
as compared to a J£ r 12 i 55 For example, substitution of Fe 
with the isovalent, but larger Ru atom does not intro- 
duce any additional charge into the system and is sug- 
gested to simulate the effect of 'chemical pressure', al- 
luding to the possibility of observing a T — > cT-phase 
transition in this substitution series. For the substitu- 
tion series SrFe2- :r Ru :E As2 (0 < x < 2), a significant but 
monotonous contraction of the lattice parameter c was re- 
ported, though no clear phase transition from a T — > cT- 
phase was observed^ The lack of a phase transition in 
SrFe2- a; Ru a ;As2 could be a consequence of substitutional 
disorder in the samples and the model of a chemically in- 
duced pressure volume effect is most likely an oversimpli- 
fied assumption^ Nevertheless, we wanted to discern the 
idea of Ru substitution acting as a 'chemical pressure'. 
To this end, we have collected the pressure dependence of 
the structural parameters a, c, c/a and unit cell volume 
for a 10% Ru-doped sample SrFe1.3Ruo.2As2 (see Fig.[T]). 
We have chosen a rather small Ru content to keep the 
impurity/disorder effects to a minimum. At room tem- 
perature and pressure, SrFc1.sRuo.2As2 is in the T-phasc 
with a c/a « 3.1. Similar to the parent compound, this 
sample is also metallic and non-superconducting at am- 
bient pressure. Upon application of pressure, the sys- 
tem also shows a transition to a cT-phase, though the 
transition pressure is shifted downwards (« 9±1 GPa) 
as compared to SrFe2As2 (please refer to Table. 3 and 
4 in supplementary information for more details). A fit 
to a Birch-Murnaghan equation of state was done for the 
tetragonal and the collapsed tetragonal phases. The ob- 
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taincd equilibrium volume, bulk modulus Bq and its pres- 
sure derivative B' for the tetragonal phase (0-7.3 GPa) 
are V = 191.2 A 3 , B = 63.3 GPa and B' = 2.85, respec- 
tively, and for the collapsed tetragonal phase (9-20 GPa) 
are V = 183.1 A 3 , B Q = 83.6 GPa and B' = 2.86, respec- 
tively. Comparing the experimental results of SrFe2As2 
and SrFe1.8Ruo.2As2, we can now perceive that the 10% 
Ru substitution did indeed act as a 'chemical pressure' in 
the sample and reduced the external pressure needed for 
inducing a T — > cT-phase transition in SrFe1.sRuo.2As2 
to 9 GPa as compared to the 10 GPa pressure needed for 
the (Ru-free) SrFe2As2 sample. 



4- Calculations: SrRu2As2 

The bonding scenario for Ru in the end member of the 
substitution series, SrRu2As2 is similar to that of Fe in 
SrFe2As2- Moreover, SrRu2As2 is non-magnetic at ambi- 
ent conditions in contrast to the paramagnetic SrFe2As2- 
Analyzing the effects of external pressure in SrRu2As2 
will therefore isolate all other parameters and provide a 
purely chemical picture of the transition. The c/a ra- 
tio of SrFc2As2 is rs 3.15 at ambient conditions, and is 
ps 2.65 for the cT phase. On the contrary, SrRu2As2 
has a ratio c/a ps 2.68 at room temperature^ and am- 
bient pressure, and therefore is already in the cT-phase 
at ambient conditions. Therefore, the effect of external 
pressure on SrRu2As2 must be much weaker compared 
to SrFc2As2. Collected in Fig. [3] arc the changes in the 
lattice parameters a, c and c/a obtained from LDA cal- 
culations as a function of volume. Both a and c decrease 
rather monotonously with the decrease in volume and do 
not show any intermediate anamolous expansion of the 
a lattice parameter. Moreover, in contrast to the cal- 
culated results obtained above for SrFe2As2 and in ac- 
cordance to the generally accepted trend, LDA slightly 
underestimates and GGA (Fig. 2. in supplementary in- 
formation) slightly overestimates the equilibrium volume 
with respect to the experiment ! 12 i 57 

Since this phase transition is not connected with a 
change in the symmetry, the nature of the transition can 
be first, second order, or even a continuous crossover. 
Moreover, the transition might be first order at T = 0, 
but if the critical end point terminating the cT(P) phase 
boundary line is below 300 K, one would only observe a 
continuous crossover at 300 K. The signatures expected 
for a first order transition, a second order one, and a 
crossover are respectively a jump, a jump in the deriva- 
tive, and an S-shaped behavior in the relevant property, 
in the present case, the c/a ratio. The data presented in 
Fig. Q] do not evidence a sharp jump, instead they clearly 
suggest a jump in the derivative d(c/a) / dT, especially for 
pure SrFe2As2- Thus the present data indicate the tran- 
sition in SrFe2As2 to be of second order type, instead of 
the well-established first order type in CaFe2As2-— This 
difference between CaFe2As2 and SrFe2As2 is in line with 
LDA calculations for the difference in the evolution of the 



total energy as a function of decreasing volume and c/a, 
where a second local minimum in c/a is much more pro- 
nounced in CaFc 2 As2 than in SrFe2As2-— 



B. Chemical bonding 

Ternary compounds AT2X2 crystallizing with 
the tetragonal ThCr 2 Si2 type of structure^ are 
numerou o 59 ' 60 and have been the focus of experimental 
and theoretical studies in solid state sciences since 
several decades^ This research originated with in- 
vestigations on the Si and Ge-based family of AT2X2 
compounds , 58 ' 62 ' 63 and later embraced pnictogens as X 
elements^ - — That pnictogen-based members of this 
family of compounds are always good for a surprise has 
already been recognized with the discovery of supercon- 
ducting LaRu2P2 (T c = 4.1 K) more than two decades 
ago^ The crystal structure comprises of T2X2 layers 
with edge-sharing TA4 tctrahedra parallel to the ab 
plane. These layers are additionally separated by planes 
of metal atom A (as can be seen in Fig. [7]). Systematic 
structural and chemical observations early on have 
resulted in classifying the AT2A2 compounds adopting 
the ThCr 2 Si2 structure type into two branches: (i) One, 
with a three-dimensional network built up of tetrahedral 
TXi layers, held together by X — X bonds along the 
c axis (two apex X atoms from 2 adjacent layers form 
the bond), with the A atoms embedded between these 
layers. This structural peculiarity has been mainly 
observed for silicides and germanides^ This branch is 
usually described as the ThCr 2 Si2 or CeGa 2 Al 2 type, 
(ii) The other branch exhibits a rather two-dimensional 
layered structure with large separations between the 
X atoms along the c axis. This branch is described 
as the TlCu2Se2 type. Interestingly, depending on the 
transition metal T, compounds containing pnictogens 
(well documented for phosphides) have been found to 
belong to both branches and to intermediate cases. 
Variation and substitution of the transition metal T 
(increasing the number of d electrons) or the metal 
atom A have been discovered as a means of tuning 
this structural peculiarity by causing a shrinkage of the 
X — X distances and bringing X species in close bonding 
contacts. Furthermore, it had also been discerned that 
the geometrical constraints due to the changes in the 
sizes of the constituent atoms alone was not adequate 
to facilitate the X — X bond formation. These analyses 
early on evoked questions about the electronic struc- 
ture governing the underlying bonding situation a 67 ' 70 
and the observed physical properties! 68 ' 71 ' 72 Absence 
of pnictogen-pnictogen bonds indicated by large c/a 
rations leads to a composite-like structural arrangement 
of the [T2X2] layers held together by more or less ionic 
interactions mediated by the electropositive A cations. 
Members of the TlCu2Se2 branch could therefore be 
qualitatively described within the Zintl concept^ 3 , by 
assigning a formal oxidation number -3 to the most 



6 



electronegative element (pnictogcn) resulting in a 
polyanion \T 2 X2\ 2 ~ counterbalanced by a cation A 2+ 
(e.g. Sr 2+ Fe 2+ Fe 2+ P 3 ~P 3_ ). In contrast, compounds 
of the TlCu2Se2 family exhibiting pnictogen-pnictogen 
bonds are preferentially formed by the late transition 
metals. The corresponding qualitative view in light 
of the Zintl concept is now based on diatomic [A2] 4 ~ 

entities (e.g. Ca 2+ Ni 1+ Ni 1+ [P 2 ] 4_ as discussed in 

Refs.[69l andTnl). However, regarding the T — T bonding 
interactions a slightly different trend was inferred 
from atomic distances and occupation of simple and 
idealized molecular orbital schemes, suggesting a more 
pronounced covalency for Fe-Fe bonds than for Ni-Ni 
bonds (Refs.[liS)- In a similar way of reasoning 
the chemical bonding was qualitatively rationalized in 
As-based compounds. 68 ' 71 Hoffmann highlighted this 
genuine occurrence of making and breaking of a diatomic 
(X — X) bond in the solid state for the AT2X2 struc- 
tures, and pioneered the investigation of the chemical 
bonding situation in phosphides based on extended 
Hiickel calculations, discovering structural flexibility due 
to a subtle interplay between pnictogen-pnictogen and 
pnictogcn-mctal bonding with the additional ingredient 
of packing requirements of the large A atoms 1^ 

1. Chemical bonding in k-space 

In order to investigate the development of an inter- 
layer chemical bond between As p z orbitals in SrFe 2 As 2 
as we go from the ambient pressure T-phase to the cT- 
phase, we analyze the orbital character of the band struc- 
ture ("fat bands"). Usually, fat bands are obtained by 
suitable projection of the extended wave function onto 
test orbitals with a certain character, which in a local- 
orbital code naturally are chosen to be the local orbitals 
themselves. Atom-centered projectors do not probe the 
phase relations between different sites, though. In or- 
der to obtain this phase information, which contains 
the chemical interpretation in terms of bonding char- 
acteristics, projectors are needed, which contain sev- 
eral sites. We choose projectors, which include As Ap z 
orbitals of a pair of As atoms sitting on top of each 
other across the Sr interlayer spacer. The two combina- 
tions are bonding Q a = -4= ($ x — $ 2 ) an d antibonding 

$> a * = -j= ($! + $2) within the pair. Here, the neg- 
ative sign gives bonding because of the odd parity of 
the p-orbitals. The corresponding fat bands probe the 
phase correlation between the two As p z orbitals but do 
not probe the phase between pairs of As atoms within 
the layer. Due to the unitarity of the transformation 
from orbitals to projectors the resulting bonding and 
anti-bonding orbital weights sum up to the "standard" 
orbital-projected weights, which opens not only the possi- 
bility to analyse the band structure but also the projected 
density of states according to bonding characteristics. 

Fig. 3] shows the a and a* fat bands for the ambient 
pressure T-phase and for the high pressure cT-phase of 



SrFe2As2- In both phases the bonding part of the 4p z 
bands (green or light grey) are mostly occupied. In the T- 
phase the antibonding bands are significantly occupied, 
while they get pushed up close to and above the Fermi 
level in the cT-phase. There is one antibonding band 
remaining below the Fermi level, however the total weight 
of occupied antibonding bands is strongly reduced for the 
cT-phase as compared to the T-phasc. 

In order to quantify the bond strength we plot the dif- 
ference of the integrated a and a* projected density of 
states (DOS) along with the projected DOS in Fig. [5] 
In an idealized case of completely occupied bonding and 
completely unoccupied antibonding states this difference 
of energy-resolved occupations will be monotonously in- 
creasing up to a value of two (spin degeneracy) at the 
chemical potential, signifying that two electrons form the 
bond, and will monotonously decrease with increasing en- 
ergy away from the Fermi level, to reach zero if all bands 
arc filled with a total number of four electrons. The value 
of this difference occupation at the Fermi level is an es- 
timate of chemical binding in terms of excess bonding 
electrons (EBE). Fig. [5] shows that the number of EBE 
is only 0.4 in the T-phase but is reaching nearly 1 in the 
collapsed phase. It turns out that the number of EBE in- 
creases by « 0.7 immediately before the collapse occurs, 
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FIG. 4. (Color online) Bonding (a, green or light gray) and 
antibonding (a* red or dark gray) As 4p z "fat bands" for 
SrFe2As2. The Fermi level is at OeV. Upper panel: T-phase. 
Lower panel: cT-phase. 
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FIG. 5. (Color online) The a (dashed, green color or light 
grey) and a* (solid, red color or dark grey) density of states 
for the T phase (upper panel) and for the cT-phase (lower 
panel) of SrFe2As2. The solid (black) line shows the number 
of excess bonding electrons, EBE (see text for explanation). 
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FIG. 6. (Color online) Fermi surfaces for the tetragonal T (up- 
per panel) and collapsed tetragonal cT (lower panel) phase. 
The color code shows the magnitude of the Fermi velocity. 



showing that the driving force of the collapse is the onset 
of Ap z interlayer bond formation. 

Returning to the band structure, besides the increased 
band width in the cT-phase compared to T-phase, there 
are additional changes occuring at the Fermi level. In 
the T-phase there is a strongly antibonding band start- 



ing at -0.5 eV at the M-point going to 3 eV at the T-point. 
This band is hybridized with the Fe 3d x 2_ y 2 orbitals. In 
the cT-phase this particular band now begins at around 
leV above the Fermi level and has a larger dispersion go- 
ing up to 7eV. Most importantly, this band that formed 
a Fermi surface in the T-phase, gets removed from the 
Fermi level for the cT-phase, which results in a change 
of the Fermi surface topology (Fig. [5]). In general the Fc 
bands do not change as drastically as the As bands across 
the transition. Most notably the Fe band width is not 
following the strong increase of the As band width. From 
a chemical point of view, the Fe layer stays inactive in the 
collapse. However, from a physical point of view, the Fe 
bands, especially the Fermi surface, changes from a 2D 
to a 3D system. Of the two cylindrical sheets around the 
M-point (T-phase) only one remains after the collapse, 
while of the three cylinders around T one vanishes, the 
second becomes a very small pocket, and the third forms 
a large pocket around the Z-point. This large pocket 
shows a quite 3-dimensional distribution of Fermi veloci- 
ties. The plasma frequencies change from fl xy = 2.81 eV 
and fl z = 0.96 eV in the T-phase to ft xy = 2.56 eV and 
fi 2 = 3.81 eV in the cT-phase. In fact the ratio A between 
in-plane and out-of-plane plasma frequencies goes from 
A ps 3 to A ps 0.7. Interestingly, the individual bands 
behave rather differently. The still more cylindrical band 
around M has a ratio A ps 1.2 while the pocket around Z 
has A ps 0.48, which makes this band actually tend to be 
more one-dimensional along the z axis, which is another 
aspect in the formation of interlayer As- As bonds. 



2. Chemical bonding in real space 

The electron localizability indicator (ELI, T) was eval- 
uated in the ELI-D representation according to Refs.l42l 
and [zl with an ELI-D module. The ELI-D distribution 
for the ambient condition T-phase shown in the upper 
panel of Fig. [7J has four distinct features. The valence 
(5th) shell of the Sr atoms is absent suggesting the for- 
mation of the Sr cation and the transfer of these electrons 
to the [Fe2As2] anion. The penultimate (fourth) shell 
of Sr is not specially structure d 79 ' 80 indicating that the 
electrons of this shell do not participate in the bonding 
interactions in the valence region. Between the As atoms, 
two distinct maxima of ELI-D are observed illustrating 
the absence of As- As bonds and the non-bonding (lone- 
pair-like) interaction between the neighboring [FC2AS2] 
anions. The structuring of the penultimate (third) shell 
of the Fe atoms towards the closest As is the fingerprint 
of the Fe-As bonding within this anion. For the high- 
pressure cT-phase (lower panel of Fig. [7]) in addition to 
the Coulomb interaction between the Sr cations and the 
[Fe2As2] anions, there is also bonding by the electrons of 
the penultimate (fourth) Sr shell (cf. structuring of this 
shell in Fig. |7J lower panel). This observation is similar 
to the one for Eu in EuRt^Gas^. The distinct maxima 
found close to the Fe-As contacts (ELI-D isosurface with 



8 




FIG. 7. (Color online) ELI-D of the ambient pressure T- 
phase (top panel) and the cT-phase at high pressure (bottom 
panel) of SrFe2As2- The slices in the left part reveal the 
distribution of ELI-D between the As atoms. The absence of 
As- As bonding in the T-phase is illustrated by the isosurfaces 
of T=1.27, while the formation of As-As bonds in the cT- 
phase is illustrated by the isosurfaces of T=1.18 . 

T=1.24) shows the formation of covalent bonds within 
the [Fe2As2] anion. The formation of the As-As bonds 
between the neighboring anions is visualized in real space 
by a concentration of ELI-D close to the middle point 
of the As-As contact (ELI-D isosurface with T=1.18), 
which is very similar in topology to the ELI-D distribu- 
tion in the simple example of the F2 molecule! 82 ' 83 

IV. SUMMARY 



der pressure with the lattice parameter c decreasing con- 
tinuously while the lattice parameter a increases for a 
certain pressure range. Our observation is akin to pre- 
vious reports^— on other 122 systems AT2AS2 (A = 
Ca, Ba, Eu; T = Fe, Co). From our experiments, we 
note that Ru substitution of the Fe-site works as 'chem- 
ical pressure' thereby reducing the amount of external 
pressure needed to obtain the T — > cT-phasc transition 
in SrFe1.3Ruo.2As2. Band structure calculations repro- 
duce the isostructural phase transition observed in ex- 
periments, including the anisotropic changes in the unit 
cell dimensions. Detailed analyses in fc-space (Fermi sur- 
faces, band structure and density of states) and as well 
as in real space (ELI-D) of the bonding scenario within 
and between the FeAs layers in both the T and cT-phascs 
presently provides a comprchcndiblc picture of the driv- 
ing force behind the observed lattice collapse. The ac- 
tual phase transition can be rationalized in the following 
way. Since the covalently bonded rather rigid transition 
metal-pnictogen layers are separated by the large cations 
(e.g. Sr 2+ ), short As-As bonds along the c axis cannot 
be formed due to geometric reasons in the ambient pres- 
sure T-phase. The predominantly ionic interaction be- 
tween the cations and the polyanion makes these layered 
compounds soft along the stacking c axis. Thus, apply- 
ing pressure leads to a decrease of the distance between 
the layers and therefore the As-As separation. When a 
critical As-As distance is reached to enable a sufficiently 
high enough orbital overlap of the 4p z orbitals, promptly 
the bonding interactions dominate and the phase transi- 
tion occurs. Similar to the case previously studied for 
phosphorus-based compounds ) 25 ' 32 electronic structure 
calculations outline a scenario where a stabilization of the 
As-As bonding states stabilize Fe-As antibonding states 
which have lifted up the former ones to energies near 
the Fermi level. Consequently, the Fe-As bonds become 
weaker and the lattice parameter a increases whereas c 
decreases at the phase transition. 



In conclusion, we have studied the transition under 
high pressure from a tetragonal (T) to a collapsed tetrag- 
onal (cT) phase in SrFe 2 As 2 and SrFe1.sRuo.2As2, using 
diamond anvil cells and powder XRD measurements up 
to 20 GPa at room temperature. We observe an isostruc- 
tural phase transition from a T-phase to a cT-phase at 
10 GPa for SrFe 2 As 2 and at 6 GPa for SrFei. 8 Ru .2As 2 . 
Both materials show an anomalous compression behav- 
ior (anisotropic changes in the unit cell dimensions) un- 
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